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Abstract 

Disorder on the string theory landscape may significantly affect dynamics 
of eternal inflation leading to the possibility for some vacua on the landscape to 
become dynamically preferable over others. We systematically study effects of 
a generic disorder on the landscape starting by identifying a sector with built- 
in disorder - a set of de Sitter vacua corresponding to compactifications of 
the Type IIB string theory on Calabi-Yau manifolds with a number of warped 
Klebanov-Strassler throats attached randomly to the bulk part of the Calabi- 
Yau. Further, we derive continuum limit of the vacuum dynamics equations 
on the landscape. Using methods of dynamical renormalization group we de- 
termine the late time behavior of the probability distribution for an observer 
to measure a given value of the cosmological constant. We find the diffusion 
of the probability distribution to significantly slow down in sectors of the land- 
scape where the number of nearest neighboring vacua for any given vacuum 
is small. We discuss relation of this slow-down with phenomenon of Anderson 
localization in disordered media. 
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1 Introduction and summary 

As we know, Universe is expanding with acceleration in the present epoch of its evo- 
lution [1]. It looks like that 1) acceleration becomes noticeable at redshifts z ~ 1, and 
that 2) the magnitude of acceleration itself does not strongly depend on the redshift 
z and is extremely small in natural units. Indeed, the matter density associated with 
this acceleration is of the order 10~^^g/cm^, 120 orders of magnitude smaller than 
the Planckian density. It remains unclear what is the meaning of this huge gap in the 
spectrum of the underlying fundamental theory. 

The physical nature of the source of this acceleration is also unclear. In general 
relativity we have learned that the expansion of the universe filled with any known 
type of matter should decelerate. Does Nature try to tell us that gravity itself should 
be modified at z ~ 1 [2] or is there simply an additional degree of freedom [3] with 
a very long relaxation time and with a special equation of state eA ~ —Pa? This 
problem, arguably one of the most important one in the modern theoretical cosmology, 
is known as the cosmological constant problem. 
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An attempt to resolve the cosmological constant problem stems from using an- 
thropic arguments [4]. The logic behind anthropic reasoning is simple: universes (or 
Hubble patches) with negative cosmological constant, very rapidly collapse reaching 
an AdS (anti de Sitter) Big Crunch singularity. Universes with large positive cosmo- 
logical constant, on the other hand, expand so rapidly that the large scale structure 
does not have enough time to form. In other words, observers able to measure the 
value of cosmological constants do not appear in universes with large positive and 
negative cosmological constants, and the bare fact of our existence automatically im- 
plies that the value of the cosmological constant should be small in our Universe. The 
statement that "the cosmological constant is very small" is true since if it were wrong 
there would be no observer to measure the cosmological constant. 

We have to introduce a Multiverse (see, e.g., [5]) - an immensely large collection of 
Hubble patches with different values of cosmological constants. Following anthropic 
reasoning, the issue of the smallness of the cosmological constant of our Hubble patch 
is thus resolved. 

An important argument in favor of the anthropic principle is that such a system 
of Hubble patches can be realized on the string theory landscape of metastable vacua 
[6] populated by eternal inflation. Even then, one might be left unhappy with an 
explanation of the cosmological constant's smallness based on anthropic principle, 
since many questions remain unresolved. For example, it is still unclear how to 
properly define a gauge invariant probability for a given observer to measure a given 
value of the cosmological constant (recent discussion of the measure problem can be 
found in [7,81) or how to determine the distribution function of vacua within the 
string theory land.cape0 

Instead of following the temptation of anthropic reasoning, one could try to es- 
tablish a dynamical vacuum selection principle on the landscape such that dynamics 
of the theory itself could automatically guarantee that vacua with low positive values 
of the cosmological constant are preferable. 

One such dynamical selection principle on the landscape is based on the analogue 
of the Anderson localization [10-12]. In condensed matter theory, Anderson localiza- 
tion was first discovered [13] as a phenomenon of electron wave function localization 
inside a semiconductor, provided that the disorder (for example, impurities or de- 
fects) of the effective potential that electrons are experiencing is sufficiently strong. 
Physics of Anderson localization in a semiconductor is related to existence of impu- 
rities, called localization centers, with potential so strong that they bind the Bloch 
wave function of the electron propagating in a disordered medium. 

Since Anderson localization is a general phenomenon taking place in disordered 
media, it is plausible that with a significant amount of "disorder" on the string theory 

^The latter problem as well as the problem of calculating correlation functions weighted over 
such distribution hmction turns out to be NP-hard [9]. 
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landscape may localize the wave function of the Universe in vacua with low cosmo- 
logical constant [10,11], thus providing an important dynamical selection principle. 
In this paper, we qualitatively study the effects of disorder on the string theory land- 
scape. 

Instead of solving the Wheeler-De Witt equation with disordered potential [10,11] 
and in order to capture effects of eternal inflation important for populating vacua on 
the landscape we consider vacuum dynamics equations [14] 

dP 

^ = E (r..^. - r.,^.) = Yl H^,P, . (1.1) 

J 3 

Equation f 1 1.1 1) describes eternal inflation and dynamics of tunneling between de Sitter 
vacua on the landscape. In this expression Tij are tunneling rates (inverse character- 
istic times of transition) between vacua i and j. Due to tunneling between different 
vacua on the landscape, probabilities Pi for an observer to measure a given value 
of the cosmological constant Aj in her Hubble patch {i.e., to find herself in a given 
vacuum i on the landscape) evolve as time passes. If the tunneling rates Vij have sig- 
nificant amounts of disorder, an analogue of Anderson localization may appear: there 
may exist vacua, "localization centers", such that the probability for an observer to 
live in such a vacuum is higher. 

The long time (large but finite r) behavior of the probability distribution function 
-Pj(r) is physically the most relevant to understand. We find it in three steps: 

• first we construct continuum limit of the vacuum dynamics equations (11. ip . 
This is possible when the number of vacua on the landscape is sufficiently large 
and tunneling rates T^j do not strongly fluctuate as functions of the position on 
the landscape (labeled by the index i). 

• then we average over disorder. The reason for us to do this averaging is the 
following. When disorder is weak, one can treat its effects by means of pertur- 
bation theory (see Appendix . Perturbed eigenstates and eigenvalues of the 
operator Hij in (11.11) are given by the integrals over the moduli space of the 
landscape. Averaging over the moduli space of the landscape is equivalent to 
averaging over disorder, and the latter is technically much simpler to perform 
than the former. 

• finally, we apply dynamical renormalization group (RG) methods to find the 
long (but finite) r behavior of the distribution function Pi{T) on a landscape 
with disorder. 

The main result of our paper is that the diffusion of the probability to measure 
a given value of the cosmological constant in a given Hubble patch may indeed be 
suppressed in the parts of the landscape where the number of nearest neighbors for 
any given vacuum is small. 
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This paper is organized as follows. We begin by identifying a possible source of 
disorder on the string theory landscape in the Section [2J We argue that compactifying 
the Type IIB string theory on a special class of Calabi-Yau (CY) manifolds will 
induce disorder on the string theory landscape. We consider Calabi-Yau manifolds 
which have large numbers of both long and short Klebanov-Strassler (KS) throats. 
Disordering comes from the fact that the lengths of Klebanov-Strassler throats and the 
points of their attachment to the bulk Calabi-Yau space may be essentially random. 
As a byproduct, we show that it is possible to have eternal inflation of stochastic type 
in the KKLMMT setup [15] with multiple Klebanov-Strassler throats. 

In Section[3]we derive the continuum limit of the vacuum dynamics equations (11. II) 
by only using the fact that the number of vacua on the landscape (or a given sector of 
the landscape) is very large. In Section H] we study the continuum limit of (11.11) . which 
describes the dynamics of eternal inflation on a disordered landscape. In particular, 
we determine the late time (large but finite r) behavior of the probability distribution 
Pj(r) by using dynamical renormalization group methods for various values of — 
effective dimensionality of the given island on the landscape. Section O is devoted to 
discussion and conclusions. 

2 Eternal inflation in a scenario with multiple 
Klebanov-Strassler throats 

In this Section we apply the idea of vacuum dynamics to the infiaton field in string 
theory setup. We start with the setup as laid down in [15] where the open string 
modulus on a mobile D3-brane plays the role of the infiaton field. This class of 
models belong to the D3-anti-D3-brane infiation models. The D3/D7-brane model 
of infiation [16] is another model of infiation where one of the open string moduli 
plays the role of infiaton. There are also other models of infiation within string 
theory, where some other moduli, in particular closed string moduli, play the role 
of infiaton. Some specific models include Kahler moduli infiation [17,18], complex 
structure moduli infiation [19], racetrack models [20] and axion modulus infiation [21] 
and [22] which used the a'-correction to Kahler potential [23] to produce a metastble 
dS4 vacuum instead of anti-D3-branes. 

As mentioned before, we confine ourselves to the D3-anti-D3-brane infiation mod- 
els. We consider a generic type IIB compactification on a Calabi-Yau orientifold in 
presence of fiux. Locally at few places the geometry of the compactification is known: 
these local patches are denoted by the Klebanov-Strassler throats [24]. The D3-brane 
travel between various Klebanov-Strassler throats is equivalent to hopping between 
different de Sitter vacua within the string landscape. We are interested in showing 
how disorder gets introduced once the D3-brane starts traversing between various 
throats. This begs for an actual calculation of the tunneling rate of the infiaton field 
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between various throats. The calculation of tunneling rate of the inflaton field from 
the first principle of string theory needs proper understanding of the dynamics of the 
inflaton field in the bulk Klebanov-Strassler geometry after incorporation of the back- 
reaction. We will discuss this aspect in a future publication [25]. Here we confine 
ourselves to the more qualitative feature of the problem. 

We first show that it is possible for the infiaton field to reach the value within a 
single Klebanov-Strassler throat, where stochasticity takes over the classical motion. 
This is shown by proving the existence of a robust upper bound on the value of the 
inflaton fleld in such a throat. This helps us define a probability distribution of the 
infiaton field in over the stringy landscape of vacua. In the next step, we estimate 
the maximum possible number of larger Klebanov-Strassler throats for Calabi-Yau 
compactification in presence of fiux. 

2.1 An upper bound on the inflaton fleld in a single-throat 
scenario 

In [26], it has been shown that the open string modulus which plays the role of infiaton, 
has a maximum range in warped D-brane infiation models. We will summarize their 
result here. In particular, we note that the upper bound on the open string modulus 
is independent of the details of the compact manifold and hence it is pretty much 
robust. 

Let us assume a warped compactification to four spacetime dimensions in Type 
IIB string theory setup [27], with the ten-dimensional line element 



where ?7^j, = diag(— 1, 1, 1, 1) and Qmn is the metric on a six-dimensional Calabi-Yau 
manifold Aie- Our convention for the indices are M, = 0, . . . , 9, /i, z/ = 0, . . . , 3 and 
m, n = 4, . . . , 9. The radial coordinate r of /Aq is defined by = {y^y -|- . . . -|- {y^y. 

The internal space is a Klebanov-Strassler throat or deformed conifold which is 
a cone over a five-dimensional Sasaki- Einstein base manifold Q5 [24]. This conical 
throat is assumed to be smoothly joined into the rest of the compact six-dimensional 
Calabi-Yau manifold SV^. The full metric over this compact Calabi-Yau is not known, 
however, locally at the throat it is known completely, i.e., the warp factor e^^ is known 
and it is a complicated function of r, closed string coupling Qg and fluxes. The origin 
of the coordinates r inside the throat is at the deformed tip of the conifold. We have 




,n 



(2.1) 
(2.2) 
(2.3) 
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turned on the following fluxef 



(27r)2a' / 




M 



(27r)2a' / 




H, = -K , 



(2.4) 



where M, K are integers, A is the 5*^ at the tip of KS throat, S is a 3-cycle, Poincare 
dual of A. In (12.41) the field strengths F3 = dC2 and = dB2, where C2 and B2 
are the RR and NSNS two-forms, respectively; see, e.g., the reviews [28,29] for more 
discussion. 



Figure 1: The locations tir and r of the D3 and anti-D3-brane, respectively, inside a 
single Klebanov-Strassler throat. It also shows the distance ruv from the tip of the 
throat where the throat merges with the bulk Calabi-Yau. The origin r = is located 
at the tip of the throat. 

Suppose, that the throat merges smoothly to the rest of the Calabi-Yau at some 
distance ruv measured from the tip. The anti-D3-brane is placed very near the tip, 
^iR ~ 0, see Fig. [H When rjji -C r <C r^/y, the warped throat can be locally 
approximated as AdSs x fl^, with the warp factor. 



Throat 



r=0 




uv 




(2.5) 



^We have also turned on the RR 0-form Co- 
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where R is the radius of curvature of AdSs. Since the background in (12.11) is generated 
by fluxes of amount MK, we have the relation [30] 

Here Vol(fi5) denotes the dimensionless volume of the space ^5 with unit radius. 
Under the AdS approximation of the throat, we can think of the background as the 
near-horizon geometry of D3-branes, where A^ = MK. 

If V^y' denote the total warped volume of the compact six dimensional space, 
the 10-dimensional Newton constant, the four- dimensional Planck mass, Mp, is [31] 

yv 

where k^q = |(27r)^(a')^. The warped volume of the internal space is 



Ml = , (2.7) 



= J d'yV~9e e-^^ , (2.8) 

where ge = detgmn- 

We can formally split the volume of the internal space into the volume of the 
throat, V^hroat; P^^^ ^^^^ the Calabi-Yau, SVfl", 

Ve^ = VZo.t + SV,- ■ (2.9) 

Under AdS approximation, the throat contribution is 

rruv 

VZo.t = Vol(fi5) / drr'e-'^^^^ 
Jo 

= 2Ti^gMKa'\ly , (2.10) 

where in the last step, we have used the relation (12.61) . Note that (12.101) is independent 
of Vol(fi5). Since the bulk volume is mo del- dependent and we are looking for a fairly 
robust upper bound of the infiaton field, we use the fact > V^^^^^, so that we can 
write (12.71) as 

yw 

Ml > -f^ . (2.11) 

If r denotes the mobile D3-brane radial modulus with respect to the anti-D3-brane 
inside the throat, the canonically normalized infiaton field is = y/Tsr, where T3 = 
1/ ((27r)^a'^(yf>,) is the physical D3-brane tension. The maximum allowed value of the 
infiaton field is thus 0max = V^'^uv- Hence we can write 

2 ^ TAW • l-^'-l-^J 



Mp J Ml V,^^^,, 
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Using (12. 6p and the values of T3 and kiq in (I2.12p . we get the maximum allowed 
value of the inflaton field in four- dimensional Planck units, in a Klebanov-Strassler 
throat [26] 



^max 



' 4 



< T7^ ■ (2-13) 



Mp j MK 

Thus, for example, for MK = 100, ^^ax ^ 0.2Mp. 

Now we should find whether it is possible to have eternal inflation with 0max — 
0.2Mp within this setup. The distance r between the D3 brane and the anti-D3 brane 
after an appropriate canonical normalization = plays the role of the inflaton 

in the four-dimensional effective field theory. The effective potential for the inflaton 
has the form 

VU) -^-^(1- (2 14) 

where 0o = VTs^o denotes the location of anti-D3-brane. Around 00; ^-C-, near the 
tip of the throat, the geometry deviates from AdSs x 5*^. From (I2.14p . we get the 
slow-roll parameters 

e(</>) - Mp(}iy^^MM (2 15) 

If the conditions e, 77 <^ 1 are satisfied, then we have an inflationary dynamics in our 
setup. Among the two conditions, (12.161) is the more restrictive one. From there, we 
have 



> 



/ 20 N V6 



Therefore, if the vacuum expectation value of is large enough, then we have infla- 
tionary dynamics. 

During one Hubble time At ~ H^^ the value of the inflaton field decreases by 

A i ; A ^ dV Ml dV Mp /— — Ml 

A0~0At —- -J2e((h) = (2.18) 

Now let us take the quantum fluctuations of the inflaton field into account. During 
the same Hubble time the quantum fluctuations 50 of the field are generated; for 
the fluctuations with the wave length / ~ ~ we have the following estimation 
for the amplitude 



H_ / 2V 87r0^ _ / Stt 0g 

2^ ^ y 37rM|, ^ y 3M^MK ~ V SMK Mp 



(2.19) 
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Therefore, if reaches the value 




) 



1/5 



(2.20) 



or higher, the stochastic force becomes more important than the classical and the 
inflation enters the eternal regime [32]. Comparing fl2.13p with (12.201) . we find that 
the maximum value of the inflaton field allowed by string theory in a single Klebanov- 
Strassler throat, is at least ten times bigger in magnitude than what is needed to 
trigger the stochastic fotcei Since the bound obtained in (12.131) is much robust and 
independent of the details of the throat, we conclude that it is possible to have eternal 
inflation even with a single Klebanov-Strassler throat. 

Finally, let us note that the regime of eternal brane inflation is in a sense opposite 
to the one usually considered in literature (see for example [34] and references therein): 
regime of deterministic slow-roll or DBI inflation corresponding to relatively small 
values of and ending with a D3-anti-D3 pair annihilation at the tip of the throat. 
Sufficiently short KS throats where D3-brane quickly approaches the tip of the throat 
effectively play a role of sinks on the multithroat sector of the landscape that we 
consider. Once the D3-brane enters such a throat, eternal inflation comes to the end 
for an observer living on the brane. 

2.2 An upper bound on the number of large 
Klebanov-Strassler throats 

In this section we will be interested in estimating the number of longer Klebanov- 
Strassler throat in a generic Type IIB Calabi-Yau flux compactification. In particular, 
we would like to get an idea about the maximum number of such longer throats. A 
statistical analysis for the number and length of such throats following the principles 
laid down in [9] was done in [35]. Our discussion will be more general than this 
analysis. 

Let us assume that the compact Calabi-Yau manifold has = 'Ylii number of 
large and number of small Klebanov-Strassler throats. The radii of curvature 

Ri of the former class of throats are large compared to the latter ones, r^. We also 
assume that Nj ^ as the larger throats have small curvatures and hence the 
classical supergravity analysis is valid. From (12.71) . (12. 9p and (I2.10p . it follows that 

^According to [33], the maximum possible value ^max of the inflaton field in brane inflation 
scenarios with warping is such that regime of eternal inflation is not allowed. The estimation of 
'/'max by X. Chen et al. is based on the rough estimation of the volume of the Calabi-Yau as 
Vjj™ ^ and subsequent requirement that R < L (as follows from the value of the 4-dimensional 
Planckian mass). The estimation of the CYq volume (j2.10|) . ()2.1H) taking warping of the throat into 
account is more accurate. 
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Larger Throats (Average size R) 




Smaller Throats (Average size r ) 



Figure 2: A Calabi-Yau manifold with many larger throats and a few smaller throats, 
the average radius of curvature of the larger throats is R. On an average they merge 
with the bulk Calabi-Yau at a distance fuv, measured from their tip. The average 
radius of curvature of the smaller throats is f. 



Ml 



Vol(fi5 
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f2.211 



Here (ruv)i and (ruv)a are the distances of the larger and smaller throats, respectively, 
measured from the tip of the respective throats, where they merge smoothly with the 
rest of the Calabi-Yau manifold. From (12.211) . it follows that 



^X:N.fiJ("uv)?<A4. 



(2.22) 



For estimation, let us assume that all the larger throats have almost same radii of 
curvature, i?j ~ ^ (see Fig. [2]). Moreover, we also assume that their merging points 
in the bulk of the Calabi-Yau are almost the same and equal to some average value 
(ruv)j — Tuy. Since = Nj is the total number of larger throats, we obtain 



N; 



UV ^ ^^-'p 



(2.23) 



Using (12. 6p and the relation 0max = V^'^UV) where 0max is the average value of the 
maximum bound of the inflaton field inside the larger throats, we have 



where M, K are some average values of the fluxes inside the larger throats. 

The larger throats have smaller curvature, so in their presence we can trust the 
semi-classical supergravity analysis. Their presence makes the dynamics of the in- 
fiaton over the string theory landscape more interesting. In particular, one major 
dynamical problem is to find out the tunneling rates of the infiaton field between var- 
ious larger Klebanov-Strassler throats. One can address this problem in the language 
of four-dimensional effective field theory. In general, the bulk Calabi-Yau geometry 
can be quite complicated and we expect that the classical traversing of the D3-brane 
from one throat to another may be impossible. In that case, we can look for the 
tunneling rate due to barrier penetration following references [36-38]. 

3 Continuum limit of vacuum dynamics equations. 
Tunneling graph 

In the previous Section we have identified one source of disorder on the string theory 
landscape: a sector on the landscape corresponding to a certain class of manifolds 
string theory can be compactified to, namely, Calabi-Yau manifolds with a number 
of warped Klebanov-Strassler throats attached to the bulk part of the Calabi-Yau 
manifold (see Fig. [2]). Since KS throats can be attached randomly to the bulk CY 
and have random lengths, dynamics of eternal inflation on the sector of the landscape 
will be affected by this randomness. 

Let us now consider the vacuum dynamics equations (11.11) on a generic landscape 
(or a sector of the landscape). 

To simplify the discussion, consider first a part of the landscape which consists 
only of de Sitter vacua, i.e., vacua with positive cosmological constants. We want to 
understand the dynamics of tunneling between de Sitter vacua. 

The vacuum dynamics equations (II. ip define the dynamics of probabilities Pi to 
measure a given value of the cosmological constant in a given Hubble patch (z.e., for 
an observer to find herself in a de Sitter vacuum i). In general, the index i enumerates 
all possible de Sitter vacua on the landscape or just vacua within a single (sub)sector 
of the landscape we are interested in. For example, i can label vacuum states with 
different number of fluxes such as in the Bousso-Polchinski landscape [39] or different 
KS throats in the sector of the landscape corresponding to tunneling between separate 
KS throats such as in the case discussed in Section [21 

The tunneling rates between de Sitter vacua i and j are given by Vij. These 
tunneling rates are not necessarily given by the Coleman - de Luccia instantons: for 
a example, on the sector of the landscape considered in Section [2], Vij are given by 
inverse characteristic times of D3-brane travel between different KS throats. Many 
other sectors on the string theory landscape may exist due to a non-trivial dynamics 
of volume and complex structure moduli of the CY manifold. In the following, we 
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do not need to know the tunneling rates explicitly. By studying the properties of 
the solutions of vacuum dynamics equations fll.ip will let us draw lessons of physics 
interest which hold in general. 

The general solution of the system fll.ip can be represented as a sum [40] 

P.(r) = J2 , (3.1) 

i=i 

where Aij are constants of integration, nmax is the total number of vacua on the 
landscape and ujij are complex (in the general case) functions of tunneling rates 
Tij. More precisely, the system of rimax first order differential equations (11.11) can be 
transformed to a single n'^^^ order differential equation 

Ed-' Fi 
a,^ = . (3,2) 

i=i 

In (13. 2p ttj are known functions of the tunneling rates T^j. Then, the rates Uij are 
given by rimax roots of the n^^^^ order polynomial 

?^max 

^ aj(-)^tj^^^ = , V A; = 1, . . . ,nmax • (3.3) 
i=i 

Clearly, when the number of vacua rimax is very large finding the spectrum of cUj/s 
by directly solving the polynomial equation ( 13. 3p is hopeless. However, when the 
tunneling rates Vij are weakly fiuctuating functions of the position i on the landscape 
(and nmax ^ 1) the vacuum dynamics equations (II. ip have a continuum limit. In this 
limit we can adopt powerful methods of the theory of partial differential equations 
and functional analysis. 

To determine how the continuum limit of the vacuum dynamics equations (11. ip 
emerges, we notice that (II. ID define dynamics of tunneling (hopping) on a graph (see 
Fig. [3]): each vacuum is a node in the graph denoted by the (composite) index and 
links between different nodes have weights given by tunneling rates Tij. As we have 
depicted in Fig. [3], parts of the graph can be isolated. These are called islands on 
the landscape (see also [40,41]). The Hausdorff dimension of a particular island is 
related to the number of nearest neighbors in such a tunneling graph. For example, 
for vacua represented by red dots in Fig. ^ N = 1. On the other hand, N = 2 for 
the vacua represented by green dots in Fig. [31 

For simplicity, we will first consider the quasi-one-dimensional case = 1 (two 
nearest neighbors for any given vacuum on the island). In this case, exactly two 
adjacent vacua exist for any given vacuum i on the island and tunneling rates to 

^On the "multithroat" sector of the landscape discussed in the Subsection |221 it is not very large, 
but still we can take Umax of order of 100. 
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Figure 3: The tunneling graph of vacua. The island of vacua denoted by green is 
well separated from "red" vacua because the tunneling rates r25, are small. The 
dimensionality of the island of "red" vacua is 1 (two nearest neighbors exist for 
any given "red" vacuum), while N = 2 for the island of "green" vacua. 
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other vacua from i are suppressed. Following the terminology introduced in [12] such 
islands are denoted as quasi-one- dimensional islands on the landscape. The system 
of equations f 1 1.11) reduces to 

= rj+l,i-Pi+l ~ ^i,i+lPi + ^i-l,iPi-l ~ ^i,i-lPi ■ (3-4) 

First of all, let us suppose that tunneling rates Fj j+i are slightly fluctuating near 
some average value f , so that they can be represented as 

r,,+i = f + 5Tl^, + 5F^+i , (3.5) 

where 



1 

2 

is symmetric part of the tunneling rate and 



'5^m= o(^r,,,,+i + 5Fi+i,) (3.6) 



1 

^r"i+i = - (5Fi,i+i - (5Fi+i,,) (3.7) 

is its antisymmetric part. Apparently, this situation can be realized in the multi- 
throat scenario we have discussed in Section [2] - for this, it is necessary that the 
lengths of different KS throats and the distances between them in the interior CY 
were not very much different. 

Let us rewrite the system (13. 4p as 

dP 

= f (P,+i - 2P, + P,_i) + A , (3.8) 

ctr 

where 

A = 5F^+^(P,+i-P,)-5F^_i(P,-P,_i)+ (3.9) 
+<5F^+i,(P,+i + P,) - 5F^_l(P, + P,_i) . (3.10) 

It is apparent from (13.81) that the continuum limit of (13. 4p exists and it is realized as 
a diffusion equation 

BP P 

when A is set to zero. For non- vanishing A 7^ 0, (13.111) has additional terms and it 
reads 

^^^l(^iD + W'(„)f^)^2l(W'(n)P), (3.12) 

where SD"^ = 6r^6n'^, dD"" = ST'^dn'^ and all higher order terms with respect to Sn 
are neglected. 

The physical meanings of symmetric and antisymmetric fluctuating parts STfj, (5F°^- 
of the tunneling rate are clear. Equation (13. 4p describes a Brownian motion of a 
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single "particle" in a random medium|^ Symmetric part of the tunneling rate ST^j 
corresponds to a fluctuation of the diffusion coefficient D, while antisymmetric part 
5r"j to a random force acting on the particle while propagating in the medium. Since 
symmetric and antisymmetric fluctuation parts are independent we can separately 
analyze the "random force" and effects associated to the fluctuations of the diffusion 
coefficient D. 

Let us now include AdS sinks to the discussion and also take into account the 
effects coming from finite volumes of expanding de Sitter bubbles. To understand the 
dynamics of the volume- weighted probability [7,40,41] we need to add extra terms 
SHiPi (if T is the world time of the four-dimensional observer) to the right-hand side 
of (11. ip . The continuum limit of the vacuum dynamics equations for a quasi-one- 
dimensional island on the landscape with AdS sinks then takes the form 

^ = ^{(D + <5/^^(n))|^) + 2|- {6D%n)P) + {3H{n) - P(n))P , (3.13) 

where r^{n) is the rate of tunneling from a given vacuum n to an AdS sink. 

By performing similar analysis for quasi-2d (A^ = 2), quasi-3d (A^ = 3) etc. islands 
on the landscape one finds the continuum limit of the vacuum dynamics equations 

dP „ . „ A / 9 /_„9P\ d 



dr 



if AdS sinks and the volume effects of dS bubbles are neglected and 

in the general case. The main difference from the quasi-one-dimensional case (A^ = 
1) is that the diffusion coefficients 6D'^ and 6D'^ are now vectors; the diffusion of 
the "particle" is anisotropic (with 6D^ being the fluctuating anisotropic diffusion 
coefficient), while the random force 6D'^ acting on the "particle" is a vector. Effects 
due to the fluctuating diffusion coefficient and random force will be treated separately 
in the next two sections. 

As a warm up before studying effects of randomness on the landscape, let us briefly 
discuss continuum limit of Bousso-Polchinski (BP) landscape [39]. The BP landscape 
consists of a discrete set of dS vacua with different values of the cosmological constant 
A defined by the formula 

A = -A, + ^Y,qlnl , (3.16) 

k 



^Technically, this means that the motion of the particle is determined by the Langevin equa- 
tion ^ = 77(71, r) + F{n), where 77 is a random force with correlation properties of a white noise 
{r]{n,T)r}{n,T')) — 2{D + 6D''{n))S{T — r'). Equation ()3.12p is the corresponding Fokker-Planck 
equation. 
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where Aq is a bare vacuum energy, are the numbers of units of a given flux k, and 
are the associated topological charges. 

The tunneling rates between different dS vacua are written in terms of Coleman 
- De Luccia (CDL) instantons [38,42-45] 

r,, = A,,e-«- , (3.17) 

where Aij ~ 1 and Bij are corresponding CDL actions. For adjacent vacua the action 
-Bj^j-i can be estimated [45] a^ 

~ • (3.18) 

8 nfqf 

The formula (13.18^ is only valid when the tunneling exponent is large, i.e., at small 
energies, when the values of cosmological constant in adjacent vacua are relatively 
small compared to the Planckian density [45]. As we can see from (13.181) . the transi- 
tions between adjacent vacua happen faster at higher rii, while at low energies they 
proceed slower. Another important feature of the BP landscape is that the transition 
downwards (towards lower energies) is always more probable than the transition up- 
wards - from lower to higher energies, because tunneling rates for adjacent vacua are 
related as 

Ti_i., = Ti^i^i exp (Si - Si_i) = Ti^i_i exp{5Si), (3.19) 

where 

Si ~ 1/A, (3.20) 

is the Gibbons-Hawking entropy. 

Let us consider transitions between vacua different by the number of quanta of a 
single flux rij — > rij — 1. For the symmetric and antisymmetric parts of the tunneling 
rate Fjj between adjacent dS vacua we have for g^nj <^ Aj, 



= ^(r^.-i + r.-i,) ^ A ■ exp (-i?M-i) • (i + l^s,^ 



(3.21) 



and ^ 

r-i-i = - F,_i,i) ^ A exp (-5i,i_i) ■ -5Si, (3.22) 

where we took 6 Si <^ 1. 

Therefore, the continuum version of the vacuum dynamics equations for this sub- 
sector of the BP landscape (tunneling between vacuum states different only by the 
number of quanta of a single flux) is 



— = A—iA- e-^- 
or on 



^For simplicity we consider flat spacetime limit. 
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This equation describes diffusion under the action of external "force" F = Ae~^''^~^SSi 
directed from higher rii (vacua with larger values of the cosmological constant) towards 
lower rii (vacua with lower values of the cosmological constant). 

Equation (13.231) describes a single quasi-one-dimensional (A^ = 1) sector of the 
BP landscape. It generalizes straightforwardly to other possible transition channels 
between vacua with different fluxes turned on. Entire BP landscape is described by a 
similar equation with N ^ 1 and all possible transition channels taken into account. 

4 Effects of disorder on the landscape 

In this Section we analyze effects of disorder on the landscape. In particular, we study 
solutions of the vacuum dynamics equations for the non-volume-weighted measure 
with disorder built in the tunneling rates Tij. We are especially interested to find 
how this disorder affects long time (large but finite r) behavior of the probability 
distribution P{n, r) for a 4-dimensional observer to measure a given value of the 
cosmological constant A„ on the sector of the landscape discussed in the Section 
EJ To determine this behavior, we pursue the following strategy. First, we suppose 
that tunneling rates Tij are not too strongly fluctuating functions of the position 
on the landscape (index i labeling vacua on the island); since the number of vacua 
on the landscape (or an island on the landscape) is considered to be very large, we 
use continuous version of the vacuum dynamics equations fl3.14p derived in Section 
131 Second, we study asymptotic properties of the general solution of the vacuum 
dynamics equations by averaging over disorder present in tunneling rates Tij and by 
using dynamical renormalization group methods. 

4.1 Random environment 

We begin analyzing (13 . 1 21) and (I3.14p by studying influence of the "random environ- 
ment" (fluctuating antisymmetric parts of tunneling rates F" ) on the dynamics of the 
probability distribution P{n, r) in the presence of disorder on the landscape. Namely, 
we consider the A^- dimensional diffusion equation 





where the "random force" 6D'^ has the correlation properties 









(4.3) 



of the white noise. Here a is a constant denoting disorder strength. 
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Such correlation properties correspond to a multithroat sector of the landscape 
discussed in the Section |21 In particular, the absence of any disorder corresponds to 
the situation when KS throats have equal lengths and are attached equidistantly to 
the bulk CY. Weak disorder corresponds to slight fluctuations of their lengths and 
distances between them. 

Strictly speaking, condition (14.31) is not applicable, for example, for the BP land- 
scape, where transitions towards smaller n are always more probable. This gives rise 
to the correlation properties 

(5D"(n)) ^ (4.4) 

to be discussed in the future publication. 

We study properties of the solution of (14.11) using methods of dynamical renor- 
malization group analysis0 

Equation (14.11) (with correlation properties (14. 2p . (14. 31) of 5D'^) can be considered 
as a Fokker-Planck equation describing a random walk of a particle in a random 
environment (characterized by the coefficients 5D'^{n)). Indeed, if the particle is 
moving according to the Langevin equation 

^ = r/(r) + 5D«(n) , (4.5) 

where the noise f/(r) is correlated as 

{Tf'{T)ri'^{T')) = D6''"''6{r-r') (4.6) 
= 0, (4.7) 

then (14. ip defines the probability to find the particle in the position n at a given time 
r. Apart from the random noise fj the particle is affected by the force in any 
point of the medium. 

To analyze random walks in a random environment, we will closely follow [46]. 
Let us hence introduce the Martin-Siggia-Rose (MSR) generating functional [47] 

Z[J] = IVPV6D^6(^^-DAP-2J2^^{6D^,P)^ 
X exp drd^nJ^T, n)P{T, n)^ exp i~ J ^ 



N 



n- — {A.i 

a 



and average over disorder present in 5D^. 

By introducing Lagrangian multipliers and transforming the functional (5-function 
into the exponent, the latter can be rewritten in terms of MSR field variables P and 



^Non-gaussian correlation properties different from (j4.2[) give rise to irrelevant terms in the ef- 
fective action (see ()4.12p below) and do not affect late time behavior of the probability distribution 
P(n,T). 



19 



P [47,48] 

Z[J] = jvPVPV6D-exp{S^sii[P,P,SD'^]) 



X exp ( / drd^n J{t, n) P{t, n)]exp{- I d^n ^^^ ^ \ , (4.9) 



where 



S'^,,^[P,PJD^]=^ j dTd^np(^-DAP-2f2^WP) \ ■ (4.10) 



d=l 

Integrating out the Gaussian disorder SD"- we find 

Z[J] = J VPVP exp(5,"g[P, P])exp (^j drd^n J{T,n)P{T,n)^ , (4.11) 

where the effective action reads 

S^ff = I [ drd^nP - dap] (4.12) 



n, r 



-^E / dTdT'd''nP{n,T)—P{n,T)P{n,r')—P{ 
^ j-^ J orid dud 

d=l 

The bare propagator of the field P in the momentum representation is 

go(^,g)= . ! n 2 > (4-13) 
—luj + Dq^^ 

where uj and q are the canonically conjugate momentum variables of r and n, re- 
spectively. Hence in the absence of disorder, the probability distribution P{n, r) is 
spreading out according to the linear diffusion law 

(jf{T)) = 2NDt. (4.14) 

In order to find the late time behavior of P{n, r), one should, in principle, work it 
out from (14.121) perturbatively in a using standard diagrammatic techniques. How- 
ever, using simple scaling arguments we can determine the late time behavior of the 
probability distribution P{n, r). From the quadratic part of the action, by multiply- 
ing momentum qhj a factor of /, we immediately find that the effective frequency u 
rescales as u' = uol^ = uP, while the Fourier components of the MSR fields P{uj, q) 
and P{uj., q) rescale as PC, and P^, where C^C, = Z^"*"^. The same argumentation for 
the interaction term proportional to a in the effective action (I4.12p . shows that the 
disorder strength a renormalizes as a' = t^'^ a. Therefore, disorder should be irrel- 
evant for islands on the landscape with N > 2, and the linear diffusion law (I4.14p is 
not affected by disorder at long time scales. 
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(a) propagator 



(b) vertex -o 



iq 

P(q,co) P(q,-co) 

co=o 

iq 

Figure 4: The MSR propagator and the vertex proportional to the strength of the 
disorder. 

To understand what happens at iV < 2 we have to generahze these naive scahng 
arguments into complete renormalization group analysis (expansion in the disorder 
strength a near N = 2). Diagrammatic rules following from the form of the MSR 
action f l4.12p are given by the bare propagator iGo{u!,q) = {u — iDq^)~^ and the 
vertex ~ —a/2 represented in Fig. m 

At one-loop level the propagator remains unchanged [46], but the vertex renor- 
malizes as 

^ = (2 - iV)a - i-a^ + 0{a') . (4.15) 
The RG flow for the dynamical exponent z is given by 

^i.2^_L.'^,0(.^). (4.16) 

(T = 0, z = 2 is a fixed point of the renormalization group flow. A second-order 
analysis shows that it is the only fixed point for N > 2 [46]. Therefore, disorder 
does not affect the late time asymptotics of the probability distribution P{n, r) for 
large A^, and the mean-square displacement {fP{T)) behaves according to the linear 
diffusion law fl4.14p at large r. 

At A^ = 2 — e, long time evolution is determined by a new fixed point 

a = A7re + 0{e^) (4.17) 
z = 2 + 2e^ + 0{e^) . (4.18) 

The mean-square displacement 

{n\r)) = I d^e— G(.;,g-)|,.o (4.19) 

behaves as 

(n^r)) ~ r'/^ ~ r^"^' . (4.20) 
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For N < 2 the diffusion process is slow, and e = 1 (quasi-one- dimensional islands on 
the landscape) is a special case treated separately below. 

Similarly, one can show that for N = 2 the mean-square displacement behaves as 



{n'{T))=ADnT[l + r^ + ...j , (4.21) 

where Dji is the renormalized diffusion coefficient, i.e., diffusion is linear with a small 
logarithmic correction. 

For N > 2 and at weak disorder a 0, the long time behavior is determined by 
the trivial fixed point a = corresponding to the usual linear diffusion law. 

It is interesting to note that the system of islands on the landscape in fact in- 
troduces a physical meaning to the RG e-expansion. This is because the number 
of nearest neighbors for a vacuum in a given island can fluctuate within the island. 
The Hausdorff dimensionality of the island (see Section [3]) can easily have non- integer 
values if the tunneling graph of the island is fractal-like. In this case, the long time 
behavior of the mean-square displacement is given by 

for 1 < < 2, and the diffusion of the probability distribution is anomalous. 

The case = 1 turns out to be the most interesting. If we naively substituted 
e = 1 into the formula fl4.20p . we would expect the diffusion process to be frozen 
even for weak disorder. This turns out to be too premature. Indeed, Sinai showed by 
using ergodic methods [49], that the diffusion process does not completely stop but 
proceeds instead very slowly according to the logarithmic law 

(n2(r)) ~ logV . (4.23) 

Notice that this behavior is independent of the disorder strength (the result holds for 
an arbitrarily weak disorder) and the disorder potentially 

The overall conclusion we come to in this Section is the following. For islands 
on the landscape with iV < 2 of the tunneling graph, the diffusion of the probability 
to measure a given value of the cosmological constant in a given Hubble patch is 
suppressed according to the formula (14.201) for any amount of the "environment" 
disorder a. More precisely, it is suppressed logarithmically for quasi-one-dimensional 
islands (A^ = 1) on the landscape and as a slow power law for 1 < < 2. For quasi- 
two-dimensional islands (A^ = 2), the linear diffusion law acquires small logarithmic 
corrections, and for N > 2 the diffusion of the probability distribution function P 
evolves according to the linear diffusion law. 



®Non-exponential behavior of i-i(r) was also reported in [50] in the different context (there the 
situation was considered when transition probabihty between vacua depends on the age of each 
vacuum) . 
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4.2 Random anisotropic diffusion 

We now turn to the discussion of effects due to fluctuations of tlie anisotropic diffusion 
coefficient . Consider the following diffusion equation 

where fluctuating anisotropic diffusion coefficients have the correlation proper- 
tiefl 

= A5,„,,5(n - n') (4.25) 
k^D^M) = 0. (4.26) 

In principle, it is necessary to consider only the case SD"^ < D, so these correlation 
properties should be modified. However, in the renormalization group analysis this 
modification will introduce higher order {i.e., irrelevant) terms in the effective action, 
and we will neglect them in the present paper. 

The Martin-Siggia-Rose generating functional now has the form 

Z|J] = / VPVW'S - f ^ ((Z, + SDiin))l^fl (4.27) 

X exp drd^n J{t, n)P{T, n)^ exp ^— J d^n ^'^"^ 

= JvPVPVSD^eMSLsR[P,P,SD']) (4-28) 
X exp ( [ dTd^nJ{T,n)P{T,n) \ exp i — [ d^n-^^ 



A 



where 



^msr[^, P, SD'] = ^ I drd-nP (^|^ - ^ ((D + 5/^^)^) j • (4.29) 
Since the disorder is Gaussian, we can readily integrate it out to find 

Z[.J] = j VPVPexp{S'^s[P, P])exp (^j drd^n J{T,n)P{T,n)^ , (4.30) 



^We again consider only the multithroat sector on the landscape described in the Section [2] 
Correlation properties of are different for the BP landscape {{SD^{n)) ^ 0) and will be considered 
elsewhere. 
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where the effective action is 

S'^ff = i j drd^nP - DAP ] (4.31) 



2 y (9nd Sn^ due due 

a,e 

Simple scahng analysis (analogous to that of Section 14.11) shows that A = is a 
fixed point of the renormalization group flow. Rescaling the length n by a factor of 
/ and requiring that the bare propagator stays invariant, will allow us to verify that 
the dynamical exponent (the one which determines how the frequency uj' = ujI^ is 
rescaled) is again 2; = 2, and the rescaling law for Fourier components of the MSR 
fields is the same as we found in Section 14.11 Disorder strength is renormalized as 
A' = Al~^ . We therefore conclude that a weak disorder in anisotropic diffusion 
coefficients is irrelevant for any A^. One can confirm this conclusion by applying 
exact renormalization group analysis. 



4.3 Random isotropic diffusion: Hermitian case 

According to the renormalization group analysis, random fluctuations of the anisotropic 
diffusion coefficients 6D^ do not affect the character of diffusion of the probability 
distribution P{n, r) in the limit r ^ 00 for arbitrary A^. This, however, does not 
mean that the dynamics of the probability distribution on the landscape is not in- 
fluenced by disorder 6D^ at intermediate time scales. To get some idea for what we 
might expect at intermediate r, let us consider the Fokker-Planck equation 

dP 

OT 

with the following Hermitiarj^ Hamiltonian [12] 

H = DA + ddiVin) ■ dd) , (4.33) 

where V{n) is a random Gaussian-distributed potential. This is a special case of (14.241) 
with only the isotropic piece {i.e., independent of (i-directions) of the anisotropic 
diffusion coefficients 6D^ non-vanishing. 

For Hermitian H the solution of the Fokker-Planck equation fl4.32p can be ex- 
panded in terms of a complete set of orthogonal eigenfunctions of the Hamiltonian 
(I4.33p . and all corresponding eigenvalues will be realo The solution of (I4.32p has the 

i°With respect to the scalar product (*', *) = / d^n^''(n)vl/(n). 

^^Let us return from the continuum limit (|3.12p back to the discrete set of the vacuum dynamics 
equations (j3.8[) . where the matrix Hij has a general form, i.e., it is neither normal nor (what is 
even more restrictive) Hermitian. In this case, the set of eigenfunctions of the operator Hij may be 
incomplete. 
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form (below we take V <^ D, but the generalization for arbitrary V is straightforward) 



P{n, t) = J2 CfcV^fc(n)e-^'=(^-^°) , (4.34) 



fc=0 



where ipk and E^. are eigenfunctions and eigenvalues of the following Schrodinger 
equation 

Ek 



di^l^u + ( ^ - W{n) ) = (4.35) 



and the effective potential W{n) is given by 

(4.36) 

Since the effective potential W{n) in fl4.35p has a supersymmetric form, all eigenvalues 
En are positive definite. If the "volume" / d^n of the moduli space of the landscape 
is finite, the ground state ipQ (normalized to ipo{n) = 1) is trivial and has zero energy 
Eq = 0. As r — s> oo (or more accurately, at r ^ E^-^) only the first term corresponding 
to the trivial ground state survives in the expansion fl4.34p . 

Let us now average over disorder, apply the RG machinery used in previous Sec- 
tions and reproduce this asymptotic behavior. After introducing Martin-Siggia-Rose 
functional and integrating out the Gaussian disorder 

{V{n)V{n')) = i2S{n - n') , (4.37) 

where is a constant, of the random potential V{n) we find the effective action of the 
form similar to fl4.3ip . Scaling analysis again shows that disorder of V{n) is irrelevant 
for arbitrary effective number of dimensions N . The physical reason for this is now 
clear: the Hamiltonian fl4.33p has a trivial extended eigenstate with Eq = Q defining 
asymptotic behavior of the probability distribution as r ^ oo. 

Late time (large but finite r) behavior is determined by eigenstates ipkip) with 
higher > 0. Also, if the ground state ipo{n) is not normalizable (such as in the case 
of an infinite landscape of vacua), higher k eigenstates are relevant in determining 
the late time asymptotics. 

The following picture [51] emerges from the analysis of the nonlinear sigma model 
corresponding to the Martin-Siggia-Rose functional fl4.3ip . For = 1 (quasi-one- 
dimensional islands on the landscape) all eigenstates ipk{n) with A; > are localized 
near the so called localization centers n = in the sense that wave functions ipk are 
exponentially peaked at n = n^. Localization centers correspond to the lowest local 
minima of the superpotential Vr(n). 

The width of ipk is called the localization length C,k- At A^ = 1 it behaves as 

~ E,' ~ k~' , (4.38) 
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i.e., it grows at small energies and reaches oo for E = Eq = 0. For small r, the lo- 
calized eigenstates with higher energy have significant contribution to the probability 
distribution fl4.34p . thus affecting the character of diffusion. 

Behavior of the correlation function at intermediate times can be estimated as 
follows. For simplicity, let us suppose that all eigenstates are localized near the same 
localization center rij [e.g., in the case when a single very deep minimum of the 
superpotential W{n) is relevant), so that 

Mn) ^ exp - ' ^ " . (4.39) 



We therefore have 



{in-n,f) ^ ^ ^ ^ ^ , (4.40) 



JdEg{E) e^+2(E)e-^" 



(4.41) 



fdEg{E)^^{E)e-^- ' 

where Ei is the energy of the first excited stat^ and g{E) is the density of states 
and it is proportional to 

g{E) oc E— . (4.42) 

For = 1 and r <^ Ei^ we find 

((n - nif) ~ . (4.43) 

This means that the diffusion is slowed down in one dimension due to the presence of 
localized states. The slow-down is not very strong, since the density of states grows 
with decreasing energy. The density of states g{E) is depicted in Fig. [5] for arbitrary 
N. 

Similarly for N = 2 (quasi-two-dimensional case), all eigenstates ipki"^) with k > 
are localized near localization centers but the localization length now grows much 
faster a.s E —>■ 

^{E) ~ e^/^ . (4.44) 

However, the density of states now remains constant as E —>■ 0, so diffusion slows 
down more effectively. 

Finally, for N > 2 there exists a mobility edge E* such that all states with E^ < E* 
are extended, while states with Ek > E* are localized. The localization length in this 
case behaves as 

~ (4.45) 

(E - E*)^ 



^^It is inversely proportional to the "volume" of the moduli space of the island and determines 
the longest time scale in the problem. 
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(b) N=2 




g(E) (c) N>=3 




E* E 

Figure 5: Densities of states g{E) for = 1, = 2 and > 3. In the = 1 case 
(a) all states ip{E,n) are localized (gray region); discreteness of the spectrum should 
be taken into account for small E. Also, in the N = 2 case (b) all states are localized, 
and the density of states remains constant. Finally, for A^ > 3 (the case (c)) there 
exists a mobility edge E* such that all states with E > E* are localized, while states 
with E < E* are not. 



near the threshold E = E*. Linear diffusion asymptotics 

{n^ir)) ~ r (4.46) 

is also reproduced when r ^ (E*)^^ ^ E^^. The diffusion process of the probability 
distribution is less effective at {E*)^^ ^ r ^ E^^: mean-square displacement behaves 

as 

((n(r) - nif) ~ . (4.47) 
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5 Conclusions and discussion 



In this paper, we systematically study dynamics of eternal inflation on the string 
theory landscape in the presence of disorder. Due to a non-trivial dynamics of moduli 
fields, many sources of disorder may exist in different sectors of the landscape. We 
identify one such sector with built-in disorder in Section [21 This sector is established 
by compactifying Type JIB string theory on a class of Calabi-Yau manifolds with a 
number of warped Klebanov-Strassler throats attached to the bulk part of the Calabi- 
Yau (see Fig. [2]). KS throats can be attached randomly to the bulk CY and can have 
random lengths. 

Eternal inflation on such sector of the landscape has a mechanism different from 
the one realized on the Bousso-Polchinski sector [39] . Our four-dimensional Universe 
is localized on a D3-brane located in one of the KS throats, and effective value of 
the cosmological constant for the four-dimensional observer living on the brane is 
essentially determined by the distance r between the D3-brane and the tip of the KS 
throat (more accurately, this distance determines the value of the effective inflaton 
field = y/TsT detected by the four-dimensional observer). The D3-brane travels 
in the throat, as a result, four-dimensional Hubble patch, the one our observer lives 
in, inflates - this is the KKLMMT scenario [15] in a nutshell. Eternal inflation (as 
we proved in Section [2], it is possible in a KKLMMT setup) corresponds to strong 
fluctuations of the inflaton field, while the distance r of the D3-brane from the tip 
of the throat is large. The D3-brane can leave the given KS throat and enter other 
KS throats, with different lengths, i.e., different scales of the cosmological constant 
measured by a four- dimensional observer. By construction, KS throats (with random 
lengths) are randomly attached to the bulk CY, and this disorder influences dynamics 
of eternal inflation. 

After identifying a possible source of disorder on the string theory landscape, we 
consider the vacuum dynamics equations (11.11) describing dynamics of eternal inflation 
on a disordered sector of the landscape (or entire landscape). One particular goal of 
our study is to understand how disorder influences the late time (r large but finite) 
behavior of the probability Piir) to measure a given value of the cosmological constant 
Aj in a given four- dimensional Hubble patch. 

Discussing eternal inflation on the landscape, usually one is interested in behavior 
of Pj(r) at r ^ oo. If the set of eigenf unctions of the operator H is complete, 
Pj(r — > oo) behavior is determined by the eigenvalue Eq of H with lowest possible 
real part. One notes that it always exists since the spectrum of H is bounded from 
below for the compact landscape. Generally, Eq is complex, and its real part can be 
negative. If this is so, P{n, r) itself as well as correlation functions weighted over this 
probability distribution are never gauge-invariant, since they explicitly depend on r. 
If the operator H is Hermitian and supersymmetric (such as for the case discussed 
in the Section IT3|l . the ground state corresponds to Eq = 0, and the asymptotics 
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P(r — i> oo) is stationary. For the case of isotropic diffusion discussed in the Section 
14.31 the form of the lowest eigenstate is trivial (constant). At weak disorder behavior 
of P(r = oo) does not depend on initial conditions for eternal inflation on such a 
landscape, on reparametrisations of r and disorder (see the Section I^TTl) . 

Although it is physically important to know Pi(r = cxo) asymptotics on the land- 
scape, finite time behavior of Pi(r) is also of great physical relevance. Indeed, because 
the moduli space of the string theory landscape is huge albeit compact, the late time 
asymptotics P;(r = oo) is reached at time scales roughly estimated as r ^ ri ~ nmax, 
where nmax is the number of vacua on the landscape, i.e., the volume of its compact 
moduli space. This time scale is extremely long, so it is plausible that our own Hubble 
patch was generated in the epoch r ^ ri 

To determine behavior of Piir) at large but finite r in the presence of disorder on 
the landscape, our strategy was to first construct continuum limit of the vacuum dy- 
namics equations (11. ip and then determine the asymptotic behavior of the probability 
distribution Pi{T) by applying dynamical renormalization group methods. 

It is easy to understand why studying the continuum limit of the equations (11.11) 
is appropriate if we consider a landscape with Hermitian H. In this case, dynamics 
of Piir) is completely determined by real eigenvalues and the form of eigenstates of 
the operator H (see the Subsection 14. 3p . At intermediate time scales r <^ E^^ (here 
El is the first excited eigenstate) the spectrum of H can be approximately considered 
continuous since the energy gap between eigenstates is of the order of inverse volume 
l/'"'max of the moduli space of the landscape (or the part of the landscape under 
consideration) . 

The use of dynamical RG methods is appropriate once recalling why RG is appli- 
cable in a generic quantum field theory (QFT): in — > limit only renormalizable 
interactions or, in other words, only relevant terms in the effective action of the QFT 
survive. Similarly, for a time-dependent problem only relevant terms survive in the 
effective Martin-Siggia-Rose action [47] describing the dynamics of Piir) (see Section 
H- 

The results of our weak disorder analysis are the following. In presence of disorder 
on the landscape 

1. large r behavior of the probability distribution Pi('r) is strongly modified on 
islands of the landscape with < 3, where A^ is the Hausdorff dimension of the 
tunneling graph corresponding to the given island. There exists a non-trivial 
fixed point in the RG flow corresponding to finite disorder (see Subsection 14.11) 
which defines this anomalous behavior. 

^•^In this respect, let us again emphasize that we study behavior of the non-volume-weighted 
measure at r — > c». Behavior of the volume- weighted measure is very different: for example, it 
is much more likely that our Hubble patch is generated at late times if probabilities are given by 
integrating over volume- weighted measure. 
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2. large r behavior of the probabihty distribution Pi{T) for islands with > 
3 is not affected by weak disorder since the only fixed point of the RG flow 
corresponds to the absence of any disorder (again see the Subsection 14. ip . 

3. on sectors of the landscape, where H is Hermitian, finite r behavior is always 
affected by weak disorder for arbitrary A^, and diffusion of the probability Piir) 
is generally suppressed for any A^. The reason of this suppression is the existence 
of localized eigenstates of the operator H, i.e., Anderson localization on sectors 
of the landscape with disorder. 

The situation in the presence of strong disorder on the string theory landscape 
is much harder to understand qualitatively. Infinitely strong disorder corresponds to 
a fixed point of the renormalization group flow of the MSR action. The physics of 
infinite disorder case is such that the diffusion of the probability distribution function 
Pj(r) is suppressed completely: as r ^ oo, the mean-square displacement (n^(r)) 
approaches a constant for arbitrary A^. This fixed point is believed to be unstable 
with respect to inverse disorder corrections. At A^ = 1 this follows from the Sinai's 
result [49] showing that the mean-square displacement diverges as r — > cxd for large 
but finite disorder. Any proof of this statement is absent for islands with A^ > 1. 
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A Why average over disorder? 

In this Appendix we will explain why it is necessary to average over disorder to de- 
termine the late time dynamics of the probability distribution, giving us a way to 
measure a given value of the cosmological constant in a given Hubble patch. To un- 
derstand the physical meaning of averaging over disorder, let us consider a particular 
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realization of the disorder on the quasi-one-dimensional (A^ = 1) landscape. This 
system is described by the continuous version of the vacuum dynamics equations 

dP 

^ = HP, (A.l) 

where the operator H is Hermitianl^ In this case, the solution of (lA.ip can be 
represented as a sum over complete set of orthogonal eigenstates of the operator H. 

Since the moduli space of the island is compact, the distribution function P{n, r) 
reaches its equilibrium asymptotics Po{n) in finite time r ~ E^^, where Ei is the 
eigenvalue of the first excited eigenstate of the operator H. The asymptotics Po{n) 
can be solved from the equation 

^(iD + 6D^in)f-^)=0. (A.2) 



dn \ dn _ 

The solution does not depend on n nor disorder and it reads 

J ""max 

The general solution of (lA.ip has the form 

P{n,r)=Y,CkMn)e-'''^ , (A.4) 

k 

where 0^ and E^ are given by the Schrodinger equation 

+ i5D^in)<P'J = Ekct>k ■ (A.5) 

The coefficients Ck are defined by the initial conditions for the probability distribution 
P(n,r = 0). 

If the disorder on the landscape is weak, i.e., 5D^ ^ D, the behavior of the prob- 
ability distribution P{n, r) at finite time can be determined by means of perturbation 
theory. Since the moduli space of the given island is compact, the spectrum of eigen- 
values Ek is discrete. In the leading approximation {i.e., in absence of any disorder), 
the spectrum of eigenstates is given by 



i?f = , (A.6) 



Yj2 

"max 



where nmax is the maximum allowed index on the island (i.e., volume of its moduli 



space). The eigenfunctions are 



\k) = <pf\n)=CkCos(^] . (A.7) 



^^Foi N > 1 this is only possible if the coefficients SD^ = 0, while the diffusion is isotropic, i.e., 
6D^ — SD ■ I, where / — diag(l, 1, . . .); the latter condition is always realized on the quasi-one- 
dimensional island. 
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In the first subleading approximation, disorder gives rise to a spectrum shift 

6eI^^ = {k\6H\k) = I dnclcos (-^^ 6H cos (-^^ (A.8) 

and the following modification in the form of eigenstates (we suppose that the spec- 
trum of states is not degenerate; the generalization for the case with degeneracy is 
trivial) 

(1) _ ^ {k'\SH\k) (0) 

^"Pk - 2^ (0) _ (0) "Pk' > i^-yj 

k'ytk ^k' ^k 



where the interaction Hamiltonian is 

I 

dn \ ^ ' On 



5H<Pk = ^{5D\n)^<pA . (A.IO) 



We find 



5E« = -cl^ I dn5D\n)sm' [^^] (A.ll 



k f , ^ „^ , . . o I kn 



n 



max 



0,4, f dn' k k'SD-in) sin {^) sin (i^) / j^,„ 



^^^^ z,..(,._,», ^^[t^] -(A-i^) 

The main conclusion we draw is the following. As we can see from the expressions 
( lA.lip and (]A.12p . after taking small disorder on the landscape into account, the 
probability distribution P(ri, r) acquires small corrections having the form of inte- 
grals of disorder over the moduli space of the island, i.e., averaging of one particular 
realization of disorder over the overall volume of the moduli space. By ergodicity, 
this is equivalent to averaging over all possible realizations of the disorder. 



B Schwinger-Keldysh diagrammatic technique 

In the Appendices [B] and [C] we mostly follow excellent lecture notes by A. Kamenev 
[52] and the classical book by L. Kadanoff and G. Baym [53], where the reader can 
find a much deeper overview of the Schwinger-Keldysh diagrammatic methods. 

The Schwinger-Keldysh diagrammatics is used to deal with quantum fields out of 
equilibrium. The crucial difference between the equilibrium diagrammatic (Feynman) 
methods and the Schwinger-Keldysh diagrammatics is that amplitudes and partition 
functions in the latter case are calculated along a closed time contour in the complex 
plane. The contour starts from t = — oo, goes to some fixed moment of time t = Iq 
and then going back to the infinite past (see Fig. [6l)l^ 

^^It is not necessary for the contour to start from infinity. It is actually not suitable to use 
such a contour for many physical situations [54]. However, using an infinitely long contour is more 
convenient, since one keeps the preparation of the initial state well under control by adiabatically 
slowly turning on the interaction term. 
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Figure 6: The closed time contour for calculating Green functions in the Schwinger- 
Keldysh formalism. 

For simplicity, let us suppose that we have a theory with a single real scalar field 
x{x). The major consequence of calculating amplitudes along the closed time contour 
is an effective doubling of degrees of freedom. In Schwinger-Keldysh diagrammatic 
technique, instead of a single field x{^) we introduce two Keldysh fields x^i.^) and 
defined on the lower and upper sides of the closed time contour. Because 
the contour is connected, the fields are not independent. In particular, the 

correlator {X^{^)x~{^)) is non- vanishing. 

After introducing the Keldysh indices a, 6 = {+, — }, the generating functional for 
the Green functions has the form 



(B.i; 



The correlation functions of the fields are given by functional derivatives with 
respect to sources at = 0. For example, the two-point correlation function of + 
and — Keldysh fields reads 

Not all of the four possible two-point Green function^ are independent. Due to 
causality constraints (see [53]) we have an identity 

G++ix,y)+G-ix,y) = G+~ix,y) + G-+ix,y) . (B.3) 

It is therefore possible to simplify the analysis of perturbation theory by doing the 
Keldysh rotation 

Xci{x) = (B.4) 
Xq{x) = -^{x^{x) -x~{x)) (B.5) 



i6/.e., the Green functions iG++ = (x+(a;)x+(a;')>> {x^ {x)x^ {x')), iG-+ = 

(x-(x)x+(x')) and iG+- = {x+ ix)x- {x')) ■ 
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since only the Green functions {Xci{x)xq{x')) = iG^{x,x'), {Xq{x)Xci{x')) = iG^{x,x') 
and {Xc\{.x)xc\{.x')) = iG^{x,x') are non-zero. G^{x,x') and G^{x,x') are retarded 
and advanced Green functions, respectively. G^ is called the Keldysh Green func- 
tion. The Green function i{Xq{x)Xq{x')) remains zero non-perturbatively because the 
identity (IB.Sp holds in all orders of the perturbative expansion [52]. 



The field Xci{x) is usually denoted as "classical" while the field Xqi^) as quantum, 
since among saddle points of the effective action there is always a solution such that 
Xq(x) = and Xci{x) satisfies the classical equations of motion [52]. 

Let us now discuss the information carried by the Keldysh Green functions. For 
a free massive scalar field we have 

G++(A;) = (e -m^ + ie)-^ -2mn{k)6{k'^ -m'^) (B.6) 
G—{k) = -{k"^ ~ + le)-^ - 2'Kin{k)6{k'^ - m'^) (B.7) 
G^-{k) = -2m{e{k'^) +n{k))6{k'^ -m^) , (B.8) 

where n{k) is the occupation number for a given mode with momentum k and 6{k) 
is a step function. 

The simplest way to verify (1B.6P - (1B.8P is to use the WKB (Wentzel-Kramers- 
Brillouin) approximation. For example, the limit x — ^ x' of the G~^^ Green function 
(neglecting the vacuum contribution) is given by 

/li^k f rJ^k 



Recalling that {x^ {.x)x^ i.^)) = iG^^{x,x) we come to the expression ( lB.6p . 
Finally, by performing the Keldysh rotation (IB.5P yields 

G^{k) = -2m{l + 2n{k))S{k^ -m^) (B.IO) 
G^{k) = (P-m^ + zO)-^ (B.ll) 
G^{k) = - Tn" - tO)-^ . (B.12) 

We conclude that the Keldysh Green function G^ carries information about the 
distribution function n{k), while the retarded and advanced Green functions G^,G^ 
give the spectrum of particles (and are independent of the distribution function n(k)). 
This separation is only valid for systems close enough to a thermal equilibrium!^ 



"'^'^Far from equilibrium the situation is not so trivial anymore, but it is possible to show that 
imaginary parts of the Green functions still carry information about the occupation numbers in the 
corresponding modes [52]. 
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C Quasi-classical Keldysh action, Martin-Siggia- 
Rose diagrammatics 

The trivial saddle point of the generating functional (IB.ip rewritten in terms of quan- 
tum and classical Keldysh fields (1B.5I1 is determined by the equations 



6S_ 
6S 



^ xq = (C.l) 
20n[xci\Xci = , (C.2) 

where Or is a retarded operator describing the dynamics of the classical Keldysh 
field, while the fields Xci, Xq were introduced in the previous Appendix. The trivial 
saddle point is S* = with Z = 1@ In fact, the condition Z = 1 holds even if one 
considers fluctuations near this saddle point [52]. 

To consider the quasi-classical limit, fluctuations of Xq should be allowed near the 
classical trajectory. Let us keep terms only up to second order in Xq in the Keldysh 
action. The semiclassical action will then take the following form 



Ss^l = 2jj dtdt' [Xq[G-TXq + {XqO''[Xcl]Xcl + cc.)] , (C.3) 



where cc. denotes the complex conjugation. 

One simple way to treat this semiclassical theory is to use the Hubbard-Stratonovich 
transformation introducing the auxiliary stochastic field ^{t) and decoupling the 
quadratic term in the quasi-classical equation. One will find that the resulting ac- 
tion is linear with respect to Xqi ^-C-, the integration over Xq leads to the functional 
5-function and to the stochastic Langevin equation 

^''[^cii^cKt) = m , (C.4) 

where 

mm) = l[G-'f{tx). (C.5) 

Another way to deal with the semiclassical theory is to integrate out the Xq field, 
since its contribution to the action is quadratic. The result is the theory of the 
classical Keldysh field 

/oo 
dtdt' Xci(0^[Xci]G^)[G^r'(G^(^''[Xci])Xci . (C.6) 
-oo 

If non-linearities with respect to Xci are neglected, the theory is a free field theory 
but with a complicated propagator. Using first quantized version of this theory, one 



^*It is possible that the Keldysh generating functional can have other saddle points (similar to 
instanton trajectories in the equilibrium situation) but they have zero contribution to the partition 
function Z. 
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can show that dynamics of the probabihty P{x, t) to find a particle excitation in the 
point X at time t, is governed by the Fokker-Planck equation. This is not surprising 
since the variation of the classical limit of the Schwinger-Keldysh equation gives the 
Langevin equation flC.4l) . 

In fact, the opposite procedure is possible: one can start from a Langevin/Fokker- 
Planck equation (or, generally speaking, with any equation describing diffusive dy- 
namics) and construct a classical limit of the corresponding Schwinger-Keldysh di- 
agrammatic technique. This procedure was first introduced by Martin, Siggia and 
Rose [47]. 

Consider the Fokker-Planck equation 

^ = HP = DAP + 6HP . (C.7) 
or 

The corresponding "generating functional" can be written as a functional 5-function 

Z= [vP6(^- hp) , (C.8) 



V dt 

which, after introducing the auxiliary field P, acquires the form 

Z = j VPVPexp (^i j dtdnP - HP^ . (C.9) 

There exists a remarkable similarity between fIC.Qp and the Schwinger-Keldysh gen- 
erating functional describing the quasi-classical approximation of the quantum non- 
equilibrium dynamics. For example, in the Martin-Siggia-Rose diagrammatic tech- 
nique, the field P plays the same role as the quantum field Xq in the Schwinger- 
Keldysh framework. 

The diagrammatic expansion of the generating functional (1C.9P is not of much use, 
since we are not interested in correlation functions of P ot P: P is itself the probability 
measure for a random walk/diffusion process governed by the corresponding Fokker- 
Planck equation. Observables in the problem under consideration (such as the mean 
square displacement (n^(r))) are related to the physics of this random walk. However, 
the correlation functions of n can be determined from the correlation functions in the 
momentum representation generated by (]C.9|) by differentiating over momenta. For 
example, we have 

(n^(r)) = J rfcue— G(c.,g-)|,-=o , (C.IO) 

where G{uj,q) is the Fourier component of the Green function (PP). This is due 
to the correspondence between the Langevin equation describing the dynamics of 
the observable n and the Fokker-Planck equation describing the dynamics of the 
probability distribution P of the diffusion process. 
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